Lower limit on decoherence introduced by entangling two spatially-separated qubits 
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It is shown that a generalization of the fluctuation-dissipation theorem places an upper bound on 
the figure of mèrit for any quantum gate designed to entangle spatially-separated qubits. The bound 
depends solely on the spectral properties of the environment. The bound applies even to systems 
performing a quantum computation within a decoherence-free subspace, but might be optimized by 
the use of non-equilibrium squeezed states of the environment, or by using an environment system 
constructed to have response confined to certain freqüències. 



Many possible systems to realise the ultimate goal 
of quantum computation ||ï[ have been proposed in re- 
cent years. In choosing systems, one has to trade ofF 
two competing requirements: one would like to be able 
to manipulate the state of the quantum bits (qubits) 
rapidly, in order to be able to perform the computation 
as quickly as possible. On the other hand, it is impor- 
tant to avoid decoherence. Typically, qubits which in- 
teract relatively strongly with thcir surroundings can bc 
manipulated more rapidly than thosc that do not, but 
also decohere more rapidly; qubits based on excitons in 
quantum dots ^, or on nuclear spins |^, provide 
examples of the fast and slow extremes respectively. In 
many cases the important quantity is the ratio of the de- 
coherence time to the characteristic manipulation time, 
often referred to as the 'figure of mèrit'. 

In this paper we make precise the connection between 
the rate at which two spatially separated qubits can be 
entangled, and the minimum amount of decoherence that 
is necessary to accomplish this, thereby placing an up- 
per bound on the figure of mèrit. This connection is in- 
evitable because spatially separated systems cannot be 
coupled directly; instead, any coupling between them 
must be mediated by some third system (the environ- 
ment). For two àtoms in a vacuum the environment is 
the electromagnètic field; for qubits in a sòlid it might 
instead be the phonons or conduction electrons. The ca- 
pacity of the environment to change state in response 
to one qubit, and the subsequent effect of this change 
on the other qubit, is what produces the coupling and 
leads to entanglement of the qubits; yet it also leads to 
the qubits and the environment becoming entangled, and 
hence to the decoherence of the qubits. (There may also 
be additional sources of decoherence, not related to the 
coupling.) It is worth noting from the outset that in or- 
der to make this connection quantitative, it is essential to 
go beyond Markovian formulations of decoherence, such 
as that of Lindblad |^ , because they take no account of 
the history of the environment. Therefore, they cannot 
represent the 'feedback' of the environment on the sys- 



tem that is the essential feature of the communication 
between the qubits. 

Consider two systems, A and i3, that encode the quan- 
tum Information. Suppose they are not coupled directly, 
but only indirectly via their mutual coupling to an envi- 
ronment E: 



H = Èa + Hb + He + \{Hea + Heb). 



(1) 



We assume the Hamiltonian is time-independent; our ar- 
gument can be generalized to the case where A and B 
are driven by a time-dependent external (classical) field. 
Let us write the coupling Hamiltonian for the A system 
to the environment as 



Hea = Va®Oa, 



(2) 



where Va is a Hermitian system A operator and Oa is 
a system E operator, and correspondingly for system B. 
(The most general form for the coupling is a sum of dif- 
ferent terms of this type, but our discussion may be easily 
generalized to that case.) 

Suppose the initial (í = 0) density matrix for the sys- 
tem consists of a direct product of some (possibly mixed) 
state of A and B with a (possibly mixed) state of the en- 
vironment E. We assume that this state of i? is a thermal 
equilibrium distribution, so we can write 



p(0) = pAs(0)(» 



exp(-^iÍ£;) 



(3) 



where (3 — í/ksT and Ze is the environment's partition 
function. Then, treating the coupling \{Hea + Heb) 
as the perturbation and working in the interaction rep- 
resentation, the reduced density matrix pab of the AB 
system evolves as 

PAB{t) = pab{0) + Ap^it) + Apf^it) + 0(A3). (4) 

The first-order term Ap^^ is given by 



Api'i(í) 



-íA / àt' 
Iq 
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{ÒA)VA{t')m+{ÒB)VB{t')m 

(5) 

where the angle brackets (...) correspond to the equi- 
librium average in the isolated environment. We keep 



-h.c. 
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Ti ^ 1 throughout. Equation (|5|) simply corresponds to 
the evolution of the AB system in the static 'applied 
field' of the environment. This evolution could equally 
well have been absorbed into that generated by Ha and 
Hb and leads neither to entangiement between A and B 
nor to decoherence (provided that any external fields en- 
tering and Hg are sufiiciently 'classical' ||^). We will 
henceforth assume that (Oa) — {Ob) — 0, and neglect 
this 0(A) term. 

To second order, we obtain 



Ap^^\t) = y2x^\ í dt' f At 



ah 

Va{t')pAB{mb{t")Cba{t'· 

dt' / dt" 

Jt' 



t') 



Va{t")Vb{t')pAB{Q)Cab{t" - t') + h.l 



(6) 



The indices {a,6} both run over A and i3, and the en- 
vironmental correlation functions Cab{t) are defined as 

Cab{t) = (Óa(í)Ófa(O)) = {é"^'Òae-'"^'Òb). The con- 
tribution to pab in equation is much more interesting, 
as the cross terms involving A and B can entangle the 
A and B systems. However, these second-order contri- 
butions also inevitably decohere the AB system, as we 
shaU see. 

We now consider a specific example in which A and 
B are single spin-1/2 qubits. We suppose that each spin 
experiences an associated classical magnètic field: 

Èa = -7Ba • a^, Hb = -tBs • à"" . (7) 

We definc ujç, to be the single-spin-flip energy in this ex- 
ternal field: üjçj = 27IBI. Initially we shall suppose that 
the applied field is in the positive z-direction. We fur- 
ther suppose that the eíFect of the interaction with the 
environment is to produce spin-flips; i.e., we choose 



Va = à^- 



Vb^oÍ 



(8) 



Note with this choice the coupling to the environment is 
perpendicular to the applied field; this will lead to a T\- 
type relaxation of the spins. Other choices of coupling 
(including more realistic ones and those that would lead 
to T2-type relaxation) are possible, but do not change 
the essential features of the result. However, we leave the 
form of the environment itself completely general at this 
stage. Re-writing the perturbation in terms of the raising 
and lowering operators (ct^Íí) — CT+e~"^°* -|-(T_e"^°*) and 
introducing two new indices s and s' which run over + 
and -, we can perform the time integrals to obtain 



Ò\pAB{^)à''As,At,u:) 



-\à1,Ò\pAB{^)<Í>s,s·{tM 



(9) 



Here the power spectrum Jab{^) is the Fourier transform 
of the corresponding correlation function 

/oo 
díC„fc(í)e''^*, (10) 

and the quantities and (j) arise from the time integrals: 

^ss'it.uj) = Tdí' /"*dí"e[-'"°(^*'+^'*")-'"(*"-*')l; 
Ja Jo 

cj)ss'{t,Lü) = í dt' í dí"e[-'"°('**'+'*'*")-'"(*"-*')l;il) 



Continuous spectrum. We first of all treat the case of 
a continuous spectrum Jabi^j) (i.e., one in which J does 
not vary significantly over a frequency range ~ t~^). For 
times that are long in comparison with the characteristic 
inverse freqüències of the environment (but still short 
compared with the characteristic motions of the qubits, 
so that the second-order expansion in A is adequate) , the 
contributions from <f> and '0 with s = s' are vanishingly 
small and we may approximate 



■KtSiuj ~ SCi'o); 



{u - swo) 

'lJjs,-s(t,Uj) « 2'KtS{Lü — SLüq). 

The second-order result thcn simplifies to 
Apfs{t)^Ap+^ + Ap^+, 



(12) 



(13) 



where the subscripts refer to the vàlues of s, s' giving rise 
to the contributions, and 



--pmisà'sjJabM 



l'K UJ — SÍÚQ 



(14) 



This can be rewritten in the effective Lindblad form 

ApAB = tí^-i[Hcs, pAB]+^(j^tiPABLl-^{LlL^, pab}^^ , 

(15) 

where the components of the effective Hamiltonian 

Hcs = H ^ + are 



H_ 



ab 



ao 



duj Jabjuj) 

— — 7 -a_a, 

27r [Lü — ujo) 

duj 1 

T^-JabKl^)-, ■ 7' 

ZTT [üJ + LÜQ) 

duj Jabiuj) ^-Buc^b-^a 



27r (oj — uJo) 



à'ÏCT^. (16) 
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Provided the Jab 7^ (i.e. provided that the fluctua- 
tions experienced by A and B are correlated), Jïeff can 
produce entanglement between A and -B, as desired. But 
there are also four Lindblad operators L^, two for posi- 
tive freqüències (linear combinations of ct"''', correspond- 
ing to correlated spin flips in the direction of the field) 
and two for negative freqüències (linear combinations of 
(t"''', correlated spin flips opposite to the field). The com- 
binations are those that diagonalize the quadratic forms 



ab 
ab 

e-'^^^A^^ J,b(^o)a^<T^ (17) 



ab 



li the environmonts of A and B are identical, symme- 
try requires that the appropriate combinations are pro- 
portional to ± . One can now see clearly the 
main point of this paper: in contrast to the Markovian 
approximation the effective Hamiltonian (generating 
the coherent evolution of the two qubits) and the Lind- 
blad operators (determining the incoherent decay) are 
not independent, because they are determined by the 
same correlation functions Job(^) of the environment. It 
is therefore in general impossible to have one without 
the other. The ratio between the two, and hence the 
maximum obtainable figure of mèrit for the gate, is inde- 
pendent of the coupling strength A to the environment; 
instead, it is determined by the spectral characteristics 
of the environmental fluctuations. 

These results are intimately connected to the 
fluctuation-dissipation theorem . The environment de- 
velops a 'response' to the applied 'field' from the first 
qubit, and it is this response that thcn drivcs the second 
qubit into a state entangled with the first. The change 
in the expectation value {Ob) when a classical perturba- 
tion Oa is applied to the environment is determined the 
frequency-dependent susceptibility 



Xba{^) = 



àuj' 



,(1 



-'^"')Jba(w') 



(18) 



Note that, since the applied field is classical, x contains 
both a positive-frequency and a negative- frequency part, 
the latter being suppressed by a factor . The real 
and imaginary parts of x connected by dispersion 
relations, and hence the dissipative (imaginary) part of 
X is directly connected to the fluctuation spectrum by 

'^{xab{í^))+'^{xba{>^)) = [JABÍuj)n{uj)+JBA{-í^)ni-uj)]. 

(19) 

Here n(w) = [exp(/3cj) — 1] ^ is the Bose occupation num- 
ber at frequency u. 

In our case, we have to distinguish between positive- 
frequency (environment transiently absorbs energy from 
System A or B) and negative-frequency (environment 



gives out energy) processes. Both contribute to the ef- 
fective Hamiltonian ([Ïg]), but each corresponds to a dif- 
ferent type of decoherence. There are therefore two sep- 
arate relations that link the effective Hamiltonian when 
the System is operated at some given frequency lu to the 
positive- and negative-frequency incoherent terms at all 
other freqüències u': 



27r Lü' — Lü 



se {+,-}• 

(20) 



A bound on the figure of mèrit for any particular gate 
can now be obtained by comparing the time needed for 
the gate operation to be generated by H^s with the de- 
coherence times governed by the Lindblad operators. 

Discrete spectrum. The opposite limit is that in 
which the power spectra of the correlation functions vary 
much more rapidly with frequency than do the functions 
Tpss'it^Lü) and (j)ss'{t,íjj). In particular, if the frequency 
response of the environment is strictly confined to dis- 
crete freqüències Í2„, we can write 



Jabii^) = E Jab.nS{üJ - íl„), 



(21) 



but must kecp the full spectral form of the functions 



(j)s,^s{t,LL!) = 
(f)s^s{t,Uj) = 
1ps,-s{t) = 



{üJ - SUJq)^ 

-t 



sin^[(cL' — suJo)t] 



1 



-1 



UJ — SUJQ {lü — SLUq)'^ 

1 



{uj — SUjQ)t 



-2isujQt 



-i{suJo+ul)t-l 



2uJo \üj — sujQ 

1 



U) + SLÜÇ) 

4 

{lo - SLUo)"^ 

-is r 1 



i{siüQ-\-iü)t 



sm 



{lü — SLÜQ)t 



1 



2ujQ \uj + SÜJQ üJ — SUJf) 

Evaluating equation (p|) for this case, we find 



ab s,s 
^b 



a°pAs(0)<T,",V.,.'(í,í^n) 

-<! <7,",(7,W(0)</'.,s'(í,í^n) +h.C. 



(23) 



This also can be written in the form of equation (|15|), 
albeit with a time-dependent effective Hamiltonian 



iíiín 



^EEE 

n ab ss' 



Jab.n 

27r 



(24) 
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and time-dependent efFective Lindblad operators: 

n ah ss' 

Note we have used the result ipss' = 'Pss' +4'*-s' -s- Once 
again the Hamiltonian and the Lindblads depend on the 
same correlation functions, although now two possible 
strategies are apparcnt that may suppress the relativa 
magnitude of the Lindblads. One is to ensure that the 
environment contains no frequency Í2„ nearly resonant 
with the transition frequency üjq] the dominant terms in 
the Lindblads then scale as í^sinc^(Aa;í), where Aui is 
the frequency interval between the operating frequency 
uq and the nearest environment frequency Í2„. A second 
strategy is to ensure that the gate's operating timescale 
coincides with the zeros of ipss'it,Lü), so that there is no 
decoherence foUowing one complete operation cycle (al- 
though there may be some at intermediate times). This 
is only likely to be possible if a single frequency Í2„ dom- 
inates the response. Note that this second strategy guar- 
antees a zero contribution to the decoherence from fre- 
quency Í7„, regardless of the magnitude of the fluctua- 
tions Jab,n- This is how a 'warm' ion vibrational mode is 
exploited to produce entanglement in an ion trap ^ ; 
our argument shows that this idea is much more general. 

Our results may be considered a generalization of the 
'fluctuation-dissipation-entanglement theorem' proposed 
by Sidles et al. |ïï|; these authors, however, considered 
just the interaction of a single spin with a harmònic bath. 
The present results confirm their expectation that such 
connections could be found for arbitrary environments. 
We also note that particular cases of the fundamental 
connection between coherent and incoherent evolution 
have been noticed before, including the Korringa law con- 
necting the relaxation rate to the Knight shift for NMR 
in metals jï^ , single-qubit manipulat ions in quantum òp- 
tics 0], and limits placed by spontaneous emission on 
quantum manipulations in ion traps |13 . 



Unfortunately, the elegant device of the Decoherence 
Free Subspace (DFS) flj is of no use in avoiding the 
type of irreducible decoherence we have described. Using 
a DFS, we would work entirely with states |'0) of the two 
qubits for which 



0. 



(26) 



However, there would then be no way to entangie the 
two qubits by the coherent part of the evolution, since 
-ffofflV') = 0. Therefore, the DFS is only useful so long as 
the dominant source of decoherence is dijferent from the 
environmental response producing the inter-qubit cou- 
pling; the DFS can reduce the decoherence from other 
sources, but not climinate the irreducible component de- 
scribed here. 

Our results suggest several ways to minimize the un- 
avoidable decoherence. Two, described above, involve 
engineering the frequency spectrum of the environment 
or carefuUy choosing the operation time of the gates. 
Another route might be by preparing special initial 
'squeezed' states of the environment, in which the fluctu- 
ations in those observables that couple to the qubits are 
reduced below their equilibrium vàlues (even, in princi- 
ple, below their zero-temperature vàlues), at the expense 
of larger fluctuations in other (conjugate) variables that 
do not couple to the qubits. Application of these ideas 
to bound the figures of mèrit for various specific types 
of inter-qubit coupling will be described in a separate 
publication. However, it should also be noted that the 
present results come from a weak-coupling expansion in 
the interaction between the qubits and the environment; 
it is entirely possible that very strong interactions in cer- 
tain Systems might lead to a different balance between 
entanglement and decoherence. 
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